Doping and Field-Induced Insulator-Metal Transitions in Half-Doped Manganites 
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We argue that many properties of the half-doped manganites may be understood in terms of a 
new two-(eg electron)-fluid description, which is energetically favorable at intermediate Jahn- Teller 
(JT) coupling. This emerges from a competition between canting of the core spins of Mn promoting 
mobile carriers and polaronic trapping of carriers by JT defects, in the presence of CE, orbital and 
charge order. We show that this explains several features of the doping and magnetic field induced 
insulator-metal transitions, as the particle-hole asymmetry and the smallness of the transition fields. 
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"Half-doped" manganites such as Rei-a^Aj^MnOs with 
X = 1/2 where Re is a 3-f rare-earth ion and A a 2+ alka- 
line earth ion have been the object of extensive studies for 
many years 0. The lowest temperature phase seems to 
be either the CE phase, consisting of ferromagnetic zig- 
zag chains with relative antiferromagnetic (AF) order (as 
in Lai/2Cai/2Mn03, where it was first proposed J3iy|, 



and in Ndi/2Sri/2Mn03 [J| or Ndi/2Cai/2Mn03 |5|) or 
an A-type phase, i.e, ferromagnetic planes with relative 
AF ahgnnient (as in Pri/2Sri/2Mn03 Q). The competi- 
tion between the CE and A phases appears even in a sim- 
ple one-orbital model because of the interplay of fer- 
romagnetic double-exchange and AF superexchange be- 
tween the core t2g spins of Mn (see also Fig. QJ. The pres- 
ence of charge and orbital order as proposed by Goode- 
nough 0] is more difhcult to establish. X-ray diffraction 
experiments do suggest the presence of large Jahn- Teller 
(JT) distortions |j,l7| with two inequivalent Mn sites. In 
the CE phase, the alternating [Zx^ — r^)/{3y^ — r^) or- 
bital order (consistent with the observed distortions) was 
shown to optimize the anisotropic hopping energy of the 
eg electrons in a more realistic two eg orbital model 0. 
The origin of charge-order was attributed to on-site or 
intersite Coulomb interactions though the latter 

tends to favor a Wigner crystal 9\ rather than the charge 
stacked order found experimentally P| . The role of the 
JT coup ling has been investigated using imposed JT dis- 
tortions as well as by extensive classical Monte-Carlo 
simulations that lead to the observed charge stacked or- 
dering [l^ . 

However, several fundamental issues remain to be un- 
derstood. One of them is the striking asymmetry with 
respect to the addition of electrons or holes. Experi- 
mentally, added electrons typically favor ferromagnetic 
metallic phases while added holes favor insulating phases 
0. In contrast, band structure arguments (8|], and treat- 
ments including JT distortions adiabatically and classi- 
cally J^j lead to metallic phases on both sides. Another 
puzzling feature, first seen in (Nd,Sm)i/2Sri/2Mn03 
and later seen to be ubiquitous , is that magnetic-fields 



~ 10—40 Tesla, which are extremely small compared with 
Neel or charge ordering temperatures ~ 200 K, induce an 
insulator- metal transition. This can be viewed as another 
manifestation of the colossal magneto-resistance (CMR) 
in doped manganites 0. An explanation is that this 
arises from the proximity of the CE phase to a ferromag- 
netic phase illJi, |14j ; but it is difficult to understand 
why the parameters in so many systems should all be so 
finely tuned as to be near the phase boundary. 

Recently, starting from a large JT coupling picture, 
a two-fluid Bg electron model, one polaronic and local- 
ized, and the other band-like and mobile, was proposed 
and shown to explain, in particular, the CMR in the or- 
bital liquid regime In this letter, we show how the 
two types of electrons can emerge from a realistic mi- 
croscopic model, even at intermediate JT couplings, in 
the half-doped case where orbital and charge order have 
to be explicitly included. Basically, they arise from a 
competition between canting of the Mn core spins pro- 
moting mobile carriers, and the JT coupling promot- 
ing polaronic, localized carriers. We show that our pic- 
ture leads to natural explanations for the particle-hole 
asymmetry around half-doping as well as the magnetic- 
field-induced insulator-metal transition at half-doping 
mentioned above. Interestingly, a similar two-carrier- 
type hypothesis was proposed based on phenomenolog- 
ical grounds in Ref. [la to understand resistivity data in 
Lai_a;Caa;Mn03 {x ~ 1/2); for which our theory provides 
a microscopic basis. JT distortions were recently tracked 
as function of field in Lai/2Cai/2Mn03 and shown to pla.y 
a crucial role near the field-induced transition 0, llSj : 
our picture is completely consistent with this. We be- 
lieve that the ideas presented here may be relevant to 
other classes of systems such as CsCeo, in which a simi- 
lar two-electron phase has been proposed ^3 ■ 

Our theory is based on the following microscopic two- 
orbital Hamiltonian for the manganites: 
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Here cj^^, creates an electron in the eg orbital a (= 
— y^, ~ ''^) in the unit cell j and a sublattice 
site labelled by a. (We use a 8-sublattice decomposition 
to accommodate the CE phase.) There are N sites and 
cN electrons with c = (1 — a;) close to 1/2. Due to a 
large Hund's coupling Jh the electron spin is assumed to 
be locked parallel to the S — 3/2 t2g core spins of Mn, 
modelled as classical vectors Sia- The hopping parame- 
ters (with 4t/3 being the hopping between (3z^ — r'^) or- 
bitals in the z-direction) include the standard Anderson- 
Hasegawa dependence on 8^, Sjb that takes care of this 
large Jh projection ■ The core spins are directly cou- 
pled by an AF superexchange, JafS'^ ~ Q.lt H is 
the external magnetic field. The last two terms include 
the vibrational energy of JT phonons (where K is the lat- 
tice stiffness of a simplified non-cooperative model) and 
their coupling to the Cg electrons. Qia and O^a repre- 
sent the amplitude and the angle of the two {Q2^ Qs) JT 
modes, and the r matrix the symmetry of their coupling 
|l2l |. On-site Coulomb interactions can be ignored in a 
first approximation when large JT distortions are present 
(as the JT coupling suppresses double occupancy) and for 
large Jh- 

We have determined the ground state of 1^ exactly 
numerically, but in the subspace of spin and distortion 
variables restricted to be periodic with a unit cell of at 
most 8 sites. This accommodates the CE state as well 
as several other competing commensurate states. Com- 
pared to earlier numerical approaches (12i | that were lim- 
ited to small clusters, our calculations are practically in 
the thermodynamic limit. We confirm the phase diagram 
that was previously obtained [12. 14J and obtain detailed 
predictions on the strength of the JT-distortions Q, etc. 
[20|. The phase diagram is given in Fig. ^and the phases 
are described in the figure caption. The strong-coupling 
phases, all insulating and charge ordered, can be under- 
stood by starting from localized Wannier orbitals cen- 
tered on alternate JT distorted sites which are fully oc- 
cupied. By virtual double exchange involving neighbor- 
ing empty sites with aligned core spins "l5| , the electrons 
gain energy in a way that depends upon the orientation 
of the JT distortion or occupied orbital 20] . A compar- 
ison of the energies of the various phases leads to the 
sequence of first-order transitions at couplings given by 
JafS^ = AtK/{9g^) and JafS^ = 8tK/{9g^) (dotted 
dashed hues in Fig. P). 

Phases that are inhomogeneous or incommensurate 
[ilj l can not be captured by the above analysis because 
of the limited size of the maximal unit-cell. We tackle 
this problem, albeit to a limited extent, by studying the 
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FIG. 1: Phase diagram of the 3D two-orbital model {T = 0, 
X = 0.5, K/t = 10). FM (resp. FMd): ferromagnetic metallic 
phase with no distortions (resp. small uniform distortions). 
FI-CO (resp. FM-CO): charge-ordered ferromagnetic insu- 
lating (resp. metallic) phase with distortions that favor oc- 
cupancy of the — orbitals. A^: ferromagnetic planes 
AF aligned with uniform distortions. A-CO: A with charge 
order. CE-CO: Ferromagnetic zig-zag chains AF ordered, or- 



bital ordered (3a; 



on the bridge sites), and 



charge-ordered (g/t > 0). G-CO: Neel AF phase with charge- 
order. Inc.: incommensurate state that interpolates between 
CE and G. Dotted dashed lines come from analytical expres- 
sions derived in the strong-coupling limit. Solid (dashed) lines 
show first-order (second-order) phase transitions. 



instabilities of the homogeneous insulating phases dis- 
cussed above with respect to particle or hole excitations 
accompanied by single site defects in their JT distortion 
pattern. For this, we find the electronic eigenvalues of 
in the presence of such defects numerically (with N up to 
1728), and calculate the energy cost or gain from filling 
the energy levels with cN electrons. 

To start with, consider the FI-CO phase at strong- 
couphng, with the electrons localized at the JT distorted 
sites with distortion Q. If we now promote a particle 
across the charge gap, it is energetically favorable for the 
JT distortion at the hole site (from which the electron is 
removed) to relax to Q — Qd- The loss in electronic en- 
ergy due to the scattering of the other electrons from the 
defect is overcompensated by the gain in elastic energy. 
The hole gets polaronically trapped, while the electron is 
mobile. Such mixed excitations thus have energies lower 
than the energies of particle-hole excitations due to band 
structure alone. For the FI-CO phase, this is demon- 
strated in Fig. [5] where, in addition to the minimum at 
Qd = (corresponding to the homogeneous phase), there 
is another minimum at Qd ~ Q, corresponding to the re- 
moval of the JT distortion at one site. Furthermore, this 
minimum softens when g/t is reduced below g/t 6.8 
(Fig. 121) , although the other homogeneous phases of Fig. 
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FIG. 2: Energy change when a single JT defect is introduced 
in the FI-CO phase. Q — Qd is the JT distortion on a defect 
site; all the other occupied sites having the same distortion 
Q. The softening of the excitation with Qd ~ Q at g/t ~ 6.8 
signals a phase transition with proliferation of defects. Finite- 
size effects are small and shown for g/t = 6.7. 



FIG. 3: Phase diagram, g/t vs. doping, x {JafS^ /t = 0.15). 
CFM: canted CE state with distortions and charge-order 
(metallic). CE trapped: CE phase with extra carriers trapped 
in JT distortions (insulating). FM: ferromagnetic metallic 
phase with no distortions. The upper curve is valid for x 
close enough to 0.5. C is a critical point ending a first-order 
line between two canted states (with different canting angles). 



^ are at higher energy at this g/t. The instability ap- 
proach therefore suggests that there might be another 
phase where such defects are energetically favorable and 
proliferate. With a small number of the above type of 
defects, a small fraction of electrons are converted from 
localized to mobile states leading to a metal with a small 
concentration of mobile electrons. This is reminiscent of 
the two-fluid picture 0|, but now extended to accom- 
modate orbital and charge order . 

We next address similar instability issues in the con- 
text of the CE phase. First consider what happens when 
the CE phase is doped with carriers. As discussed above, 
experimentally there is a strong asymmetry between hole 
and electron doping. According to de Gennes's original 
argument [i^ , canted phases are expected for small dop- 
ing (irrespective of their sign). It is known that the en- 
ergy of the fully ferromagnetic state crosses that of the 
CE state when extra electrons are added but inter- 
mediate canted phases have not been considered. They 
would naively lead to second-order transitions rather 
than first-order. We have studied such canted phases 
and obtained the optimal canting angle as function of 
X close to 1/2. Similarly to the discussion above in the 
context of the FI-CO phase, these homogeneous phases 
are in competition with inhomogeneous phases where the 
added carriers are self-trapped by JT distortions. In fact, 
we find that it is favorable to trap the added carriers at 
small dc= 1/2 — x. This leads to an insulating un-canted 
CE phase (noted CE trapped in Fig. On the elec- 
tron-doped side {x < 1/2, Fig. increasing Sc leads 
eventually to a canted metallic phase (CFM) via a first- 
order transition. This is because the JT energy gain due 
to trapping is linear in Sc, 5Etr = —Ej'^\5c\ (with Ej'^ 
obtained by solving, for all g/t, the one-defect problem 



mentioned above with one added carrier) [20|, whereas 
the energy gain from canting is quadratic, SEca — (<5c)^ 
[2^ . The latter loses for small Sc but wins for larger Sc. 
On the hole side, however, we find that canted phases are 
never energetically favorable. The asymmetry arises be- 
cause of the nature of the CE ordering. Canting leads 
to a 2-d dispersion, with a large density of states at 
the bottom of the conduction band (for electron dop- 
ing), whereas it gives a 3-d dispersion, with a vanishing 
DOS at the top of the valence band (for hole doping). 
Therefore, canting angles can indeed get large when elec- 
trons are added and compete effectively against electron 
trapping. But when holes are added, canting angles are 
much smaller and the holes get trapped by JT distortions 
for g/t > A. Hence the system remains insulating. Thus, 
our approach leads to an explanation for the asymmetry 
between particle and hole doping seen experimentally. It 
also helps us to understand why incommensurate charge 
ordered CE typ e phases seem to be favoured on the hole 
doped side HEEJ- 

An external magnetic-field applied to the CE phase 
also promotes canting. Experimentally, as discussed ear- 
lier, a field-induced insulator-metal transition occurs at 
extremely small fields. To locate the transition in our 
theory, we minimize and compare the energies of vari- 
ous 8-sublattice structures in a field, including the JT- 
distorted canted CE state, the undistorted canted state 
with the optimal (high) canting angle, etc. We find that 
for y/t < 5 the ground state switches in a first-order 
transition from a distorted canted CE phase (with the 
canting hardly changing the JT distortions) to an undis- 
torted highly-canted (or FM) phase with increasing field. 
At the transition, the system becomes metallic, there is 
a jump in the magnetization (Fig. and an abrupt 
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relaxation of all the JT distortions to zero. The transi- 
tion fields have very little to do with the magnetic energy 
scales, but are determined by the JT energies and depend 
strongly on g/t as is clear from Fig. 0] For g/t > 6.8, 
the ferromagnetic state is insulating and no insulator- 
metal transition can be found, which puts a bound on 
the values of g/t that are appropriate. In the range 
5.0 ^ g/t ^ 6.8, which may be relevant for mangan- 
ites (we need g/t > 5.0 to explain the existence of the 
A-CO phase Q| [see Fig. [T]), we find an instability of the 
distorted canted CE phase towards creation of defects, 
which suggests that the field-induced metallic phase in 
this case has the above mentioned two types of electrons. 
In all cases, there are abrupt changes in JT distortions at 
the transition, converting it from second-order (for pro- 
gressive cant ing) to first-order, in agreement with recent 
experiments [Tn. Iisj . 
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FIG. 4: Magnetization vs. field {g/t = 1-5, JafS'^ /t = 0.15). 
The first-order transitions to an undistorted highly canted (or 
fully FM) metallic phase is accompanied with a relaxation of 
the JT distortions. 

However, the transition fields obtained in our calcu- 
lations are too large compared to experiments. For in- 
stance, gfJ-sHc ^ O.lt (Fig. EI), gives ~ 140 T (with 
t ~ 0.2 eV). The discrepancy is connected with the over- 
estimation of the charge gap in our model. Three ef- 
fects need to be included to obtain a more realistic, re- 
duced estimate for the charge gap. First, the finiteness 
of the Hund's coupling, here taken to be infinite, which 
would allow for hopping even between sites with anti- 
aligned core spins; second, the cooperative nature of the 
JT phonons, causing sizeable distortions on the corner 
sites of the CE phase as well; and third, small second 
neighbour hopping. All of these would contribute to re- 
ducing the gap. Indeed, if we use the experimentally 
observed charge gap in place of the charge gap obtained 
in our model and then estimate the transition field, we 
get numbers in good agreement with observations. 

In conclusion, we have provided new theoretical in- 
sights into the physics of half-doped manganites. We 
suggest the existence of and competition between canting 



(i.e., not full ferromagnetism, which could be checked by 
neutron diffraction) induced metallicity and inhomogene- 
ity arising from the trapping of carriers by JT defects. 
It explains several features of the doping-induced (e.g. 
the particle-hole asymmetry) and field-induced insulator- 
metal transition. These ideas suggest a new two-fluid 
model with localized and mobile electrons, which extends 
the work of ref. to include orbital and charge or- 
der, which when treated with more sophisticated meth- 
ods such as DMFT can yield a satisfactory and complete 
theory of doped manganites including the regime near 
half-doping. 
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